Analytical approximation for the sphere-sphere Coulomb potential 
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A simple analytical expression, which closely approximates the Coulomb potential between two 
uniformly charged spheres, is presented. This expression can be used in the optical potential semi- 
classical analyses which require that the interaction be analytic on and near the real r-axis. 
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The semiclassical analyses of the optical potential cross 
sections play an important role in understanding the scat- 
tering mechanism of light- and heavy-ions [0. The direct 
evaluation of the semiclassical amplitudes is based on the 
calculation of the integral actions between the turning 
points of the radial motion. The turning points of the 
radial motion are solution of the equation 
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where k, A, E, and U(r) are, respectively, the wave num- 
ber, the angular momentum in units of h, the center of 
mass energy, and the complex optical potential. The in- 
tegral actions are defined by 



Sij(X) 



k(r, A) dr, 



(2) 



where and Tj are appropriate solutions of Eq. ([!]). 

In the semiclassical approximations both these quan- 
tities are treated as analytical functions of A around the 
real A-axis. This is true only if the complex optical po- 
tential is an analytical function of r, on and near the real 
r-axis. 

In the phenomenological analyses the optical potential 
between heavy-ions is written in the form 



U(r) = V(r) +iW(r) + V c (r), 



(3) 



where V(r) and W(r) are the real and the imaginary part 
of the nuclear interaction, and Vc(r) is the Coulomb in- 
teraction potential. 

As a rule, analytic functions are used for the real and 
the imaginary part of the interaction. On the contrary, 
the Coulomb part of the interaction is described using 
either the Coulomb potential of one point charged parti- 
cle with a uniformly charged sphere of radius Rc (in the 
least recent works) or the Coulomb potential between two 
uniformly charged spheres of radii, say, R\ and i?2 (in 
the more recent ones). These potentials are not analytic 
functions around the real r-axis. 

This fact, rigorously, excludes the possibility of apply- 
ing the semiclassical methods which assumes the analyt- 
icity of the potential. The lack of analyticity must how- 
ever be considered accidental. It arises from our prefer- 
ence to use simple expressions for the interaction, and we 



do not really think that it is connected with some phys- 
ical fact. Owing to this we expect that this, and similar 
troubles, can be cured with a little effort. 

To do this we here briefly recall some techniques used 
in the past to treat the point-sphere case, and we present 
a similar method for the sphere-sphere case. 

Both these Coulomb potentials can be written in the 
form 



voir) = /*(*), 
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where Z\ and Zi are the charges of the colliding part- 
ners, x = r J Rc, and fc(x) is an appropriate form factor 
depending on the case considered. For the two sphere 
case, obviously, Rc is the sum of R\ and i?2- 

The form factor for the point-sphere Coulomb poten- 
tial is given by 
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This form factor is continuous, as is its first derivative, 
on the real x-axis, but has a discontinuity at x — 1 in the 
second derivative. 

To avoid problems connected with this discontinuity, 
Knoll and Schaeffer 0], in their semiclassical analysis, 
simply described the Coulomb interaction by using the 
point-point Coulomb potential for all the r-values. This 
choice is only justified for strongly absorptive interac- 
tions, for which the contributions from the internal re- 
gion of the interaction are unimportant. In order to be 
free from this restriction, Brink and Takigawa || , in their 
uniform semiclassical analysis, preferred to approximate 
the point-sphere form factor given by Eq. |^ with 



g T (x) 



-[1 — exp(— aix — a 2 x 2 — a 3 x 3 )]. 
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By imposing the condition that the values of this form 
factor, and of its first two derivatives, coincide with the 
corresponding quantities of the exact one at x = 0, one 
obtains a\ = 3/2, a 2 = 9/8, and a 3 = 19/16. In this 
way the analytical form factor closely approximates the 
point-sphere one, with a maximum relative error smaller 
than 3%. 

A similar analytical form factor, with the additional 
constraint of being an even function of x, was also used in 
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the past. The evenness of the potential ensures that the 
scattering functions, of the multireflection expansion of 
the Brink and Takigawa approximation, have definite re- 
flection properties in the complex A-plane. These proper- 
ties are useful in the semiclassical analysis, and to ensure 
their validity in Ref. [f|, although not there explicitly 
given, the following expression was used for the point- 
sphere form factor 



§6(3 
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f/3i( x ) = — tanh(aia; + a^x 4 ). 
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(7) 



Using the same condition used by Brink and Takigawa 
one obtains for the parameters the values a± = 3/2, and 
as = 5/8. With this parameter values, the maximum 
relative error between the approximated and the exact 
form factors is lower than 4%. 
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FIG. 1. The thick curve shows the point-sphere form factor 
and the thin solid and dashed lines show the approximations 
of Eqs. Q and (Q), respectively. The relative errors of the 
approximations are shown in the inset. 

In Fig. 1 the behavior of the exact point-sphere form 
factor is compared with that of the two above approxi- 
mations. 

The use of the point-sphere Coulomb potential is 
not, however, appropriate to describe realistically the 
Coulomb interaction in the heavy-ion scattering pro- 
cesses in which the projectile and the target have similar 
radii. More than 25 years ago |-§ the role of more re- 
alistic descriptions of the Coulomb interaction between 
two heavy-ions was investigated. 

The detailed behavior of the Coulomb interaction de- 
pends, obviously, on the properties of the two charge 
distributions. In any case the potentials obtained us- 
ing realistic charge distributions can be rather well ap- 
proximated by the Coulomb potential of two uniformly 
charged spheres. This is the reason why this potential is 
used even in the more recent and detailed analyses of the 
elastic scattering between light heavy-ions [p~0|— |l3|| . 

The sphere-sphere Coulomb form factor in given by R] 



:{l-id\l-xf 
x[l - ^d(l -x)(5 + x)}} x < x < I 



x > I, 
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where b = 1 + 1/a, c = 3/5a 2 , d = (a + 1/a + 2)/4, 
xq — (a — l)/(a + 1), with a = R1/R2 and i?i > i?2- 

This formula is less popular than the equivalent ones 
given in Refs. ||[||^], but is here preferred for its simplic- 
ity. 

The sphere-sphere Coulomb form factor is continuous 
for real ir-values up to the third derivative. The first dis- 
continuities appear in the fourth derivative at x = xq and 
x=l. 

A very good approximation for the sphere-sphere form 
factor can be obtained by using a slight generalization of 
Eq. (R) 



f 



(r) = - tanh(aix + a^x* + a 5 x b ), 
x 



(9) 



where the parameters ai, 03, and a$ depends only on the 
ratio o between the larger and the smaller of the two 
radii. With respect to Eq. (^), only the fifth power term 
was added in order to preserve the evenness of the form 
factor. To give more flexibility to this form factor the 
parameters are allowed to be determined, for each value 
of a, by fitting Eq. (||) to points calculated using Eq. (||). 

This technique was used to obtain the sphere-sphere 
form factor parameters in a recent semiclassical anal- 
ysis Q of the optical potential scattering of 16 O on 
12 C at .©Lab = 132 MeV. In the optical potential con- 
sidered the Coulomb interaction was originally described 
|]l3| using the potential of two uniform charge distribu- 
tions of radii R\ — 3.54 fm and R2 — 3.17 fm. With 
these radii, the values obtained for the parameters are 
ai = 2.387, a 3 = 1.071, and a 5 — 1.683. Using our opti- 
cal code, the differences between the cross sections cal- 
culated with the exact and the approximated Coulomb 
potential are completely negligible. In fact, calculating 
the cross sections from 1° to 180° with a step of 0.25°, 
the maximum relative error is lower than 0.09% 

The need to use a fitting procedure to obtain the best 
estimate of the parameters may be considered a drawback 
of the approximation (||). This complication is however 
unnecessary because the dependence on a — 1/a of the 
best fit parameters is rather smooth and, for practical 
purposes, it is sufficient to know the values of the best fit 
parameters on a rather sparse grid of a values. 

The reliability of the fit procedure and the possibility 
of interpolating the best fit parameter values was tested 
using 101 a; values, equally spaced between and 1.5, for 
101 values of a — 1/a, equally spaced between and 1. 
In the worst case (a = 0, which corresponds to the point- 
sphere case) the maximum relative error of the best fit 
form factor is less then 0.8%. 
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FIG. 2. Values of the best fit parameters of the 
sphere-sphere approximated form factor as a function of a. 
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sphere form factor for the a values from 0.0 to 1.0, with 
a step of 0.1. The differences between the exact and the 
approximated form factors cannot be appreciated with 
the thickness used for the curves. In the inset of the 
same figure, for the a = 0.0, 0.5, and 1.0 cases, we show 
the relative errors between the exact and approximated 
form factors. 

In Tab. | we give the values of the best fit parameters 
for a values spaced by 0.1. The last column of the table 
gives the corresponding maximum relative error. Inter- 
polation methods of the tabulated values can be used 
to obtain the parameter values at intermediate a values. 
Using a linear interpolation at the middle a values the 
maximum relative error is within the bracketing error 
values for a < 5.0, and does not exceed the 0.21% for the 
other a values. 

In all the tests that we have performed, the substitu- 
tion of the sphere-sphere form factor with its analytical 
approximation given by Eq. (0), does not produce ob- 
servable effects with the thickness of the lines normally 
used to plot an optical cross section. 

In conclusion we can say that the lack of analyticity 
of the Coulomb potentials, traditionally used in the op- 
tical model analyses, can be easily cured without any 
physically appreciable effect. The same is expected to 
be true also for the, more or less, similar lacks of ana- 
lyticity in other term of the optical interactions, which 
are sometimes used. Using appropriate analytical func- 
tions it should be possible to closely approximate these 
interactions too. These analytical approximations could 
be used in the semiclassical analyses which requires that 
the interaction is expressed in terms of analytical func- 
tions. 
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FIG. 3. Sphere-sphere form factors for a from 0.0 to 1.0, 
in steps of 0.1 (curves from below to above). In the inset the 
relative errors for the a = 0.0, 0.5, and 1.0 cases are shown 
by the solid, dashed, and dotted lines, respectively. 



TABLE I. Best fit parameters and maximum relative er- 
ror for the analytical approximation ^ to the sphere-sphere 
Coulomb form factor with different a values. 



The dependence of the fit parameters on the quantity 
a is shown in Fig. ^. 

In Fig. H we shown the behavior of the exact sphere- 
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